Matching of Wilson loop eigenvalue densities in 1+1, 2+1 and 3+1 

dimensions 



Francis Bursa" 

^Institut fiir Theoretische Physik, Universitdt Regensburg, 93040 Regensburg, Germany 



Abstract 

We investigate the matching of eigenvalue densities of Wilson loops in 5*^7 (A^) lattice gauge 
theory: the eigenvalue densities in 1+1, 2+1 and 3+1 dimensions are nearly identical when the 
traces of the loops are equal. We show that the matching is present to at least second order in 
the strong-coupling expansion, and also to second order in perturbation theory. We find that 
in the continuum limit there is matching at all values of the trace for bare Wilson loops. We 
confirm numerically that there is matching in these limits and find there are small violations 
away from them. We discuss the implications for the bulk transitions and for non-analytic 
gap formation at iV = oo in 2+1 and 3+1 dimensions. 



1 Introduction 



The eigenvalues of an SU{N) matrix are pure phases, e*", and are gauge invariant for a 
closed Wilson loop. In [1] the eigenvalue distributions of Wilson loops in D = 1 + 1 and 
D = 2 + 1 SU{N) lattice gauge theories were studied, and a remarkable result was found: if 
the eigenvalue distributions of Wilson loops f/^^" of size nxn were calculated in D and D 
dimensions, with the couplings and Pd' chosen so that the expectation values of the traces 
in the fundamental representation mJJ,^" = ;^ReTrir{[/^^"} were equal, 

(«r"(/5D))D = {uT"{M)d', (1) 

then the two eigenvalue distributions would match. This was found to apply over a wide 
range of couplings, loop sizes, and values of N. Preliminary results showed that the matching 
also applied to 3+1 dimensions. Matching between eigenvalue distributions in 1+1 and 3+1 
dimensions was also found in [2j for SU{2). The eigenvalue distribution of smeared Wilson 
loops in 3+1 dimensions in the large-iV limit appears to approach the D = 1 + 1 SU{oo) 
distribution [3], suggesting the matching may also work for smeared loops, at least at = oo. 
An approximate matching of smeared Wilson loops between 1+1 and 3+1 dimensions has also 
been observed in SU (2) [4j . 

If the matching found in [1] is exact it will have implications for the 2+1 dimensional and 
3+1 dimensional theories. Firstly, in 1+1 dimensional lattice gauge theory there is the Gross- 
Witten transition [5] at = oo. This is a third-order phase transition which occurs when the 
trace of the plaquette is 0.5. At the transition a gap opens in the plaquette eigenvalue density. 
That is, on the strong-coupling side of the transition, the eigenvalue density is non-zero for all 
a, but on the weak-coupling side it is only non-zero for some range of angles — < a < ac 
with etc < vr. If the matching is exact this gap opening, and hence the Gross- Witten transition, 
would be reproduced in 2+1 and 3+1 dimensions. In 2+1 dimensions there is a bulk transition 
at this value of the trace which is however different to the Gross- Witten transition, having a 
specific heat peak which is not present in 1+1 dimensions [Ij, so if the Gross- Witten transition 
is present it would have hidden in some way by this transition. Similarly in 3+1 dimensions 
there is a strong first-order bulk transition at A^ = oo (in fact for all A^ > 5) [HI El E] , which 
would have to hide the Gross-Witten transition. 

Secondly, in 1+1 dimensional SU{oo) gauge theory in the continuum there is a non- 
analyticity, the Durhuus-Olesen transition, where a gap opens in the Wilson loop eigenvalue 
distribution at a critical Wilson loop area Acrit [HI [ID] • If the matching is exact such a non- 
analyticity would also have to be present in the D = 2 + 1 and D = 3 + 1 cases. In the 
latter case, and if it occurs at a finite physical lengthscale, it could provide an explanation 
for the rapid crossover between perturbative and non-perturbative physics that is observed in 
the strong interactions [HI [121 113] • It has also been suggested [3j that the non-analyticity at 
A^ = oo could be used to match perturbative and non-perturbative physics and thus calculate 
the string tension in terms of Aqcd- 

However, exact matching would also imply that Casimir Scaling JT^, which is exact in 
1+1 dimensions, must also be exact in 2+1 and 3+1 dimensions. Casimir Scaling has been 
found to be a good approximation for tensions between sources in various representations of 
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SU{3) in 3+1 dimensions [151 [IB], but it is known not to be exact. It is also not exact for 
fc-string tensions in both 2+1 and 3+1 dimensions [17]. Hence the matching cannot be exact, 
and must presumably be violated at some level. 

The results in [1] were not the main purpose of that work and hence were somewhat 
limited; for example, the matching of traces in eq. ([I]) was not exact and there was no error 
analysis. The purpose of the present work is to carry out a more focused, quantitative study 
of the matching and of any violations that may be present, and to assess the implications for 
gap formation. 

The structure of the paper is as follows: in the next section we show analytically that 
matching is present in both the strong-coupling and weak-coupling limits, and we consider 
the continuum limit. In Section [3] we present our numerical results. We conclude in Section |H 

2 Strong— coupling and weak— coupling limits 
2.1 Strong coupling limit 

The eigenvalue density of a Wilson loop p{a) can be expressed as a sum of Fourier components 
Fk cos{ka). The Fourier coefficients Fk are just the traces of powers of the Wilson loop: 



So if the traces of all powers of two Wilson loops are equal their eigenvalue densities will 
necessarily be equal too. 

Consider the evaluation of such traces in the strong-coupling expansion for a planar Wil- 
son loop of area A in lattice units. Apart from a possible constant term, the leading order 
contribution comes from the diagram in which we tile the Wilson loop with k sheets of pla- 
quettes, giving a contribution proportional to {up}^"^, where (up) is the trace of the plaquette. 
This is independent of the number of dimensions D, as are the group-theoretical prefactors. 
The trace in the fundamental representation is just (up)^, so the condition ([1]) is equivalent to 
matching the values of (up) . Thus to leading order in the strong-coupling expansion we expect 
that the traces of all powers of the Wilson loops, and thus their eigenvalue distributions, will 
indeed match. 

There is a possible loophole in this argument. The traces of powers of a Wilson loop can 
be reexpressed as combinations of traces of the Wilson loop in different irreducible representa- 
tions. Some of these representations (for example the adjoint) can be screened in 2+1 and 3+1 
dimensions, leading to a perimeter law instead of an area law. This should lead to violations 
of the matching. However, in practice what happens is that for small loops, the tiling with 
fewest plaquettes, which will give the leading term in the strong coupling expansion, is the one 
in which there are sheets of plaquettes tiling the plane of the loop, giving an area law. Tilings 
with sheets of plaquettes forming "tubes" around the outside of the loop, giving a perimeter 
law, only start to have fewer plaquettes for rather large loops. In the case of a square L x L 
loop, this will happen approximately when = 16(L — 1), or at about L = 15. Such large 
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loops are difficult to study numerically, and in any case will have extremely small traces: if 
we are in strong coupling, so {up) < 0.4, say, the trace in the fundamental representation of 
a 15 X 15 loop will be less than approximately 0.4^^ ~ 10~^°. Traces in other representa- 
tions will be even smaller. So the eigenvalue density will be indistinguishable from the Haar 
measure in practice, and there will trivially be matching between different dimensions. 

Higher-order terms in the strong-coupling expansion do depend on the dimensionality of 
space-time, so we expect these to cause violations of the matching. The ffist such term is 
one in which one of the sheets of plaquettes tiling the Wilson loop has a "bump" in it. This 
involves an additional four plaquettes compared to tilings in which all the sheets are flat, so is 
suppressed by a factor (up)^. However, we show in Appendix El that in fact the contributions 
from the bump are identical in all dimensions when the traces are matched as in eq. ([1]). 

So the ffist possible violations come in at the next order in the strong-coupling expansion. 
In SU{2), the next order is tilings where one sheet has a "bump" in it, all the other sheets 
are fiat, and there is an extra pair of plaquettes somewhere in one of the sheets. This is order 
(up)''^'^^. In SU{3) we must add three plaquettes instead of two, so the tiling is of order 
(Mp)*^"^"*"^, and for higher we can add another "bump", requiring an extra 4 plaquettes and 
giving tilings of order {up)^^^^. 

These contributions are very small for all but the smallest loops. The largest contributions 
are those for k = 2 and for SU{2). Assuming the prefactors to be of order one, and taking the 
largest (up) for which the strong-coupling expansion is valid to be approximately 0.4, these 
will be of order 0.001 for the plaquette, which should be easily observable, but only of order 
10~^ for the 2x2 loop, which is unlikely to be detectable, and even smaller and hence almost 
certainly undetectable for larger loops. Thus we expect that when strong coupling applies, the 
eigenvalue densities of Wilson loops will be observed to match, except for the very smallest 
loops. 



2.2 Weak coupling limit 

We now consider the eigenvalue densities of Wilson loops in perturbation theory. As discussed 
in Section [2m if the traces of all powers of two Wilson loops are equal their eigenvalue densities 
will match. Furthermore, traces of powers of a Wilson loop can be expressed as combinations 
of traces of the Wilson loop in various representations. Thus to show that two Wilson loops will 
have matching eigenvalue densities it is sufficient to show that their traces in all representations 
are the same. 

Wilson loops satisfy "Casimir Scaling" to second order in perturbation theory; that is, the 
trace of a Wilson loop in a representation R is given by 

j^TrnW = e-<'-^-^^'\ (3) 

where Cji is the quadratic Casimir for the representation R, and foi.g'^) is a representation- 
independent function (which can however depend on the dimension of space-time D). 

If eq. (II]) is satisfied it is clear the values of the functions fo and fo' must match. Then it 
follows that the traces in all other representations will match and hence the eigenvalue spectra 
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will match also. 

Thus it seems that we expect matching when the coupling constant g'^ is sufficiently small 
that second-order perturbation theory is a good approximation. However, there is a subtlety 
involving representations with numerically large Casimirs Cr. These can be important when 
g"^ is small since then the trace of the Wilson loop approaches 1 so the eigenvalue distribution 
must be narrow. Then Fourier components with large fc, that is traces of large powers of 
the Wilson loop, are important, and these can in general involve representations with large 
Cr. Now consider the trace of the Wilson loop in such a representation to leading order in 
perturbation theory: 

^'TtrW = 1 - W^CRg\ (4) 

Here Wi is a representation-independent constant which depends on the Wilson loop size and 
the dimension D. If Cr is large the quadratic term in (jlj) can also be large even though g'^ 
is small. Then we expect that higher-order terms in perturbation theory could also be large, 
including those which violate Casimir Scaling at third order and higher. 

However, consider the next term in the expansion of TirW . This is of the form: 

^TtrW = 1 - W^iC^Z + {\w^Cl - W2CR)g\ (5) 

In the limit of very large Cr the quartic piece is dominated by the ^WfC'^g'^ term. This term 
comes from the exponentiation of the tree-level term WiCRg'^: 

^-WrCns' = 1 - W.CRg' + \wlClg' + 0{g'). (6) 

The exponential will contribute a piece proportional to C|j to the g"^^ term in the perturbation 
expansion. If the remaining parts of the g^"' term are lower order in Cr, as is the case for 
the g^ and g'^ terms, the pieces coming from the exponential will dominate. The higher-order 
terms are not known, but if this is the case to all orders in (7^, which seems plausible, the trace 
of the Wilson loop will be dominated by a term of the form ([3]) and there will be Casimir 
Scaling, and hence matching of the eigenvalue densities. 

In the particular case of = cxo there is another argument for matching in perturbation 
theory which does not depend on the behaviour of representations with large Casimirs. For 
sufficiently weak coupling the eigenvalues are all near unity and then they will plausibly be 
in the same universality class as the Gaussian Unitary Ensemble, and so the density will be a 
'Wigner semicircle" . This will be true in all dimensions, again leading to matching. 

2.3 Continuum limit 

The strong-coupling analysis in Section [271] applies when the lattice coupling constant, (3 = 
is sufficiently small. This is the exact opposite of the continuum limit, for which we take /3 to 
infinity, so we do not expect the analysis to be useful in the continuum limit. Specifically, we 
expect that strong-coupling will apply below some critical coupling jSc- The trace of a Wilson 
loop of area A is approximately {up)^, where {up) oc /3 is the trace of the plaquette. Hence 
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we expect the strong-coupling analysis to apply only for traces less than which indeed 
rapidly goes to zero in the continuum limit. 



In 1+1 dimensions in the continuum Wilson loops obey an exact area law, (m^) = e'""" ^ , 
where a is the string tension, a is the lattice spacing, and L is the size of the loop in lattice 
units. So the boundary of the regime in which the weak-coupling analysis in Section 12.21 
applies corresponds to a finite value of the trace. However, this is not the case in either 2+1 
or 3+1 dimensions. In both these cases there is a "perimeter term", due to the self-energy 
of the sources which propagate around the perimeter of the loop. The leading correction to 
In(M^) is aLV{a), where V{r) is the Coulomb potential, which is proportional to g'^N\n{r) 
and g'^N/r in D = 2 + 1 and D = 3 + 1 respectively. Combining this correction with the area 
term, we get 



We see that in both cases if we take the continuum limit, keeping aL fixed, at a constant value 
of the coupling g'^N, the trace will decrease. So the boundary of the weak-coupling regimes 
will shift to lower traces in the continuum limit. This will happen only logarithmically in 
D = 2 + 1, but linearly in D = 3 + 1. 

On the weak-coupling side of this boundary the traces of the Wilson loop in different 
representations will obey Casimir Scaling. But in D = 1 + 1 there is Casimir Scaling at all 
couplings, not just in weak-coupling. And by the arguments of Section 12.21 Casimir Scaling 
implies matching, so as long as we are on the weak-coupling side of the boundary there will 
be matching with loops in 1+1 dimensions. For sufficiently large L the minimum trace for 
which the coupling is weak will become arbitrarily small, and there should be matching at all 
values of the trace. 

3 Results 

3.1 Preliminaries 

We discretise D-dimensional Euclidean space-time to a periodic cubic lattice with lattice 
spacing a and size in lattice units. We assign SU{N) matrices, Ui, to the links / of the 
lattice. We use the standard Wilson plaquette action 



where Up is the ordered product of the SU{N) link matrices around the boundary of the 
plaquette p. We simulate the theory using a conventional mixture of heat bath and over- 
relaxation steps applied to the SU (2) subgroups of the SU (N) link matrices. 

As discussed above, we expect that matching of eigenvalue spectra will only occur in 
certain limits. Specifically we expect to see matching in the strong-coupling limit, and when 
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the physical size of the loops is small enough that second-order perturbation theory is valid. 
As discussed in Section I2.3[ due to the "perimeter term" this will be the case for loops with 
any finite trace in the continuum limit in both D = 2 + 1 and D = 3 + 1. So we expect to see 
matching for any value of the trace in the continuum limit also. 

Thus we expect that, for each N, the part of parameter space in which matching does 
not occur will be constrained on three sides: in the direction of decreasing trace by strong 
coupling, in the direction of increasing trace by perturbation theory, and in the direction 
of increasing Wilson loop size (in lattice units) by the perimeter term. A fourth constraint 
is given by the loop size, which of course cannot become negative. Our strategy will be to 
determine for each Wilson loop size whether there is indeed a range of traces where matching 
does not occur, bounded by the strong-coupling regime on one side and the weak-coupling 
regime on the other, and then to increase the Wilson loop size to check whether this region 
moves to lower values of the trace, as expected due to the perimeter term. 

We begin our investigation in SU (2) since it is easiest to obtain the large statistics necessary 
to detect small violations there. We then turn to SU{3) and SU{6), to look at how things 
change with increasing and, in the latter case, to investigate possible effects due to the 
first-order bulk transition which is present inD = 3 + l [6l[7l[H]- These calculations will allow 
us to draw some conclusions about what happens at N = oo, and in particular whether the 
Gross- Witten and Durhuus-Olesen transitions will be reproduced there. 

To perform a comparison for a particular loop size n at a particular value of the trace {uw) 
we first generate configurations at couplings Prun in D = 1 + 1, D = 2 + 1, and D = 3 + 1. We 
then reweight to couplings Prew where the traces match the target value of the trace exactly. 
On every sweep we calculate Wilson loops on every 1 to 2 lattice sites, and extract their 
eigenvalues. We store these in a histogram with a large number ribin (typically 4096) of bins. 

From these histograms we estimate the first Fourier components F^, defined by F^. = 
jf Sill cos(/caj) where Aj = e*°* is the z'th eigenvalue, and the first moments fik of the 
eigenvalue distribution. The finite width of the histogram bins means that these estimates 
will have errors of C ( ) . As we shall see below, the relevant values of k are usually in 

\ bin / 

the range 2 to 20, for which this is ~ 10~^ and hence negligible. We compare the Fourier 
components and the moments of the eigenvalue distributions, searching for any statistically 
significant differences between the two distributions. We also compare the two histograms 
directly. To reduce the noise we also construct coarser histograms with ... ,4 bins, 

and compare these. 

Of course, these channels in which violations can show up are not independent; a generic 
violation will give a signal in all of them. What we find in practice is that the Fourier 
components are the most sensitive. That is, as we increase the statistics (at a loop size and 
trace for which we expect to see a violation), we first see a statistically significant difference 
in the Fourier components, and only for higher levels of statistics do the differences in the 
moments and histograms become significant. Thus from here on we will concentrate on our 
results for the Fourier components. 

Furthermore, there is also a pattern in which Fourier components show significant differ- 
ences first. We find that in SU{2), a graph of the Fourier coefficients Fk typically looks like 
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that shown in Fig. [TJ Fq is always 1 by definition, Fi is the trace and is always positive, 
and for larger k there is a rapid decrease in F^ to some maximally negative value, followed 
by an asymptotic approach to zero. We find that the Fourier components where significant 
differences first show up are those around the minimum. This is perhaps not surprising; is 
is easier to observe a small absolute difference if the quantity being observed is large. If the 
Fourier component is large and positive it is not far from 1 and thus largely perturbative. In 
this case we expect that we will not see any violations, as we discussed in Section 12. 2[ So 
indeed we expect that violations should be easiest to see for the Fourier components which 
are large and negative. 

For larger N the graph of Fourier components Fk has more maxima and minima; we then 
find significant differences first appear in the Fourier components around the first minimum. 
We will refer to the value of k at which the first minimum in F^ occurs as k'. The corresponding 
Fourier component is then F^i. Since the significant differences first appear in the Fourier 
components around the first minimum F^i is a useful quantity to use as a measure of the size 
of the violations, if any. 

3.2 SU{2) 

We now turn to our results, starting with those for SU(2). Presumably violations of eigenvalue 
spectrum matching will be largest when we are furthest from the limits analysed in Section [2] 
where we expect matching. This leads us to look at small loops at moderate values of the 
trace. We begin with the smallest loop, the plaquette. 

Our results for the plaquette are summarised in Table [H We see that, as found in [Ij, 
the eigenvalue distributions in different numbers of dimensions when the traces are matched 
are indeed very similar. We illustrate this in Fig. [2], where we plot the plaquette eigenvalue 
densities at (up) = 0.7, where the differences are largest. Despite this, we see that the 
differences are barely visible by eye. Turning to the Fourier components, which we plot in 
Fig. [3l we see that the differences are slightly more clear, but still extremely small. Note that 
these are the largest differences we see; at most values of the trace the differences would be 
impossible to see on such plots. 

We see that, as discussed at the end of Section [XT| the largest differences in Fig. E] are 
around the Fourier component at the minimum, Fk'. To investigate how the differences vary 
with the trace (up), we plot in Fig. Hlthe differences in Fk' between 1+1 and 2+1 dimensions 
and between 1+1 and 3+1 dimensions, as a function of the trace. 

Proceeding from small to large values of the trace, we see first of all that there are no 
detectable violations at (up) = 0.07, where the coupling is very strong. Violations become 
visible, possibly at (up) = 0.3, and certainly at (up) = 0.4, still within the strong-coupling 
region. After this they rapidly increase, reaching a maximum when the trace is around 0.8 for 
the difference between 1+1 and 2+1 dimensions and 0.6-0.7 for the difference between 1+1 
and 3+1 dimensions. Finally they decrease again as we approach {up) = 1. We see the same 
features for other values of k near k', and also (though with relatively larger statistical errors) 
for the moments and histograms of the distributions. 

All this is as we would expect. There is a large range of traces in strong coupling for 
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which there are no observable violations of the matching. Violations start to appear around 
the upper limit of strong coupling, due to the higher-order terms in the strong-coupling 
expansion (though they are smaller than we estimated; presumably the prefactors are small). 
Then there is an intermediate-coupling region where there is no reason to expect matching 
and where we indeed see violations. Finally the coupling becomes weak and the matching 
reappears due to Casimir Scaling. 

For the larger values of the trace, the volumes used in D = 2 + 1 and especially in D = 3 + 1 
are extremely small. Indeed, our calculations for {up) = 0.99 in 3+1 dimensions were carried 
out at a coupling (3 = 75.2, for which the lattice spacing is ~ 10^^'^ fm! Thus one might worry 
that the small differences we see there are not due to true infinite-volume differences between 
the theories but are instead due to finite-volume effects. To check for this, we have repeated 
our calculations on different sized lattices. We find no difference in the violations in matching 
as we change the lattice size, making it unlikely that these are finite-volume effects. 

We now turn to our results for 2x2 Wilson loops. These are listed in Table and the 
corresponding differences in F^i are plotted in Fig. O We can compare these to our results 
for the plaquette in Table [1] and Fig. |H Proceeding from small to large traces, we see first 
of all that there are violations for 2x2 loops at much smaller values of the trace than for 
the plaquette. This is because the boundary of strong-coupling has shifted to a much smaller 
trace — indeed at (u^) = 0.04 the trace of the plaquette is 0.447, 0.444 and 0.438 in 1+1, 
2+1 and 3+1 dimensions respectively, so we are already near the edge of strong coupling. The 
violations of matching are just becoming detectable here. They then become larger, as for 
the plaquette, but reach maxima and start decreasing at lower values of the trace than for 
the plaquette. This is presumably the effect of the perimeter term lowering the value of the 
trace at which weak-coupling begins. We see that the peak is shifted to a lower value of the 
trace in 3+1 dimensions compared to 2+1 dimensions, as expected since the perimeter term 
is stronger in the former case. 

We also see that the magnitude of the violations is smaller by a factor of 5 or so compared 
to the plaquette. This is not just because of the perimeter term, which causes larger loops to 
be at a weaker coupling at the same value of the trace. For example, comparing the plaquette 
at a trace of 0.8 to the 2x2 Wilson loop at a trace of 0.5, we see the violations are much 
larger for the plaquette, despite the fact that the couplings are roughly the same. Presumably 
what is going on is that there are relatively large lattice corrections to Casimir Scaling for the 
plaquette, which spoil the matching, but these corrections are smaller for the 2x2 loop. 

Finally we look at our results for 4x4 loops, which are listed in Table[3]and plotted in Fig. [61 
We see that the difference between 1+1 and 2+1 dimensions is non-zero for much smaller 
values of the trace than for the smaller loops, as expected. For the difference between 1+1 
and 3+1 dimensions the situation is less clear — the peak appears to be around (m^) = 0.5, 
as for 2x2 loops — but it is clear the violations are now extremely small. 

Above we have described the behaviour of the Fourier component Fk' for the 2x2 and 
4x4 loops, but as for the plaquettes we see essentially the same behaviour for other Fourier 
components with k near k', and for the moments and histograms of the eigenvalue densities. 
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3.3 SU{3) and SU{6) 

We now turn to our results for SU{3) and SU{6). Our results here are less extensive and are 
aimed mainly at establishing that the features we saw for SU (2) carry over to larger N. 

Our results for the plaquette in SU (3) are summarised in Table [Hand are plotted in Fig. [71 
We see a similar pattern to the corresponding results in SU{2): there is very good matching 
at small values of the trace, followed by violations of the matching at intermediate values, 
which decrease again as we approach (up) = 1. The violations appear to be somewhat larger 
than for SU{2). 

Our results for the plaquette in SU (6) are summarised in Tableland are plotted in Fig. [HI 
Note that due to the bulk transition there is a range of (up) around 0.5 which is inaccessible 
in 3+1 dimensions; the closest we can get is (up) = 0.435 and (up) = 0.52 on the strongly- 
coupled and weakly-coupled sides respectively. The overall pattern of violations is similar to 
that for SU{2) and S'f/(3): again we have good matching at small and large values of the 
trace and violations in between. There is no sign of anything exceptional happening around 
{up) = 0.5, the region where the bulk transitions occur. 

The maximum violations now appear quite large (~ 0.02). To judge whether this is in fact 
large, we can compare the values of Fk' at a trace of 0.5 in Tables [5] and [3 We see that they 
are around —0.04 for the plaquettes and —0.105 for the 2x2 loops, so distributions in which 
the Fk' differ by ~ 0.06 are certainly possible. So the violations observed are still significantly 
smaller than they could be. 

Our results for 2x2 loops are summarised in Tables [6] and [7] for SU{3) and SU{6) respec- 
tively. They are plotted in Figs. [9]and[T0l Again, the pattern of violations is broadly similar 
to that seen in SU{2), shown in Fig. [5] The biggest change is that the peak violations are 
larger and occur at a lower value of the trace in SU{Q). Why this should be is unclear; it 
occurs in the intermediate-coupling region where we do not have an analytic understanding 
of the source of the violations. In any case the violations are still small, and they decrease in 
both the strong-coupling and weak-coupling regions, as they should. 

The differences described above are for the Fourier component Fk', but we see the same 
patterns for nearby Fourier components and for the moments and histograms. 

In summary, we see that despite some differences in the details, the violations in SU (3) and 
SU{6) follow the same broad patterns as in SU{2), and are consistent with our expectations 
based on strong-coupling and weak-coupling analysis. 

4 Conclusions 

We have shown that the remarkable matching between eigenvalue spectra of Wilson loops 
observed in [H [2], [3l H] can be understood analytically in both the strong-coupling and weak- 
coupling limits. In the strong-coupling limit, it follows from the fact that the leading-order 
contributions come from tilings that are flat, and hence the same, in all dimensions. In 
the weak-coupling limit it follows from the fact that Wilson loops satisfy Casimir Scaling 
up to second order in perturbation theory in all dimensions. Due to the "perimeter term". 
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perturbation theory applies at all values of the trace in the continuum limit, so the matching 
will apply there also. 

However, our numerical results have shown that away from these limits the matching of 
eigenvalue densities is not exact. These differences remain small, of order 0.02 at most. Why 
this should be when we are well away from both the strong-coupling and weak-coupling 
regimes is unclear. 

Apart from their small size, the overall pattern of violations follows our expectations: it 
approaches zero in both the strong-coupling and weak-coupling limits, and has a peak at 
intermediate couplings. The trace at which this peak occurs decreases as the size of the loop 
increases, as expected due to the effect of the perimeter term. 

In the large- limit in 1+1 dimensions, the Gross- Witten transition occurs when the 
trace of the plaquette is 0.5 |5j. At this value of the trace we observe finite differences between 
the eigenvalue densities in 1+1 and 2+1 dimensions which show no sign of decreasing (indeed, 
they appear to increase) with A^. This implies that the Gross-Witten transition will not be 
reproduced in D = 2 + 1, which is consistent with the results found in [T], where a peak was 
observed in the specific heat at the bulk transition in 2+1 dimensions which is not present at 
the Gross-Witten transition. 

In 3+1 dimensions a jump in the trace of the plaquette occurs at the bulk transition 
making values around 0.5 inaccessible. However, we see differences in the eigenvalue densities 
on either side of this jump, suggesting that even if the bulk transition was removed there 
would still be significant differences at a trace of 0.5 and thus any "underlying" transition 
would again be different to the Gross-Witten transition. 

The situation is different for the Durhuus-Olesen non-analyticity in Wilson loop eigenvalue 
densities which occurs at A^ = oo in the continuum in 1+1 dimensions [9l HO]. This occurs 
when the trace of the Wilson loop is ~ 0.135. Due to the perimeter term this will be in 
the weak-coupling limit in the continuum limit, where we expect matching. Thus we expect 
that the non-analyticity will be reproduced in both 2+1 and 3+1 dimensions. However, since 
it occurs in the weak-coupling limit it is in the far ultraviolet, so it seems unlikely to be 
physically significant. 

This problem could be avoided by using smeared, instead of bare, Wilson loops, as was 
done in [21 S] . If the perimeter term was eliminated completely by smearing the Wilson loops 
would have a trace of at a constant lengthscale in physical units in the continuum limit 
(although the lengthscale would be determined by the details of the smearing procedure). 
However, the Wilson loops would then be given completely by the area term, which does not 
satisfy Casimir Scaling exactly in 2+1 or 3+1 dimensions. Then the traces of Wilson loops 
in different representations will not match and so neither will the eigenvalue densities, and so 
the non-analyticity in D = 1 + 1 will not be reproduced for smeared Wilson loops in 2+1 or 
3+1 dimensions. 

We still expect gap formation to occur, since the eigenvalue density becomes flat as ^ 
on the one hand, and forms a Wigner semicircle as -u^ — > 1 on the other. In between a gap 
must form. However, since the lengthscale it will occur at will depend on the details of the 
smearing procedure, it is unclear whether the gap formation will have any physical significance. 
Also, since the eigenvalue density will be different from that at the Durhuus-Olesen non- 
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analyticity, it is not clear what universality class the gap formation will fall into. This is a 
question that needs to be resolved in order to implement the idea of matching perturbative 
and non-perturbative physics through the transition, as set out in [3]. 

In summary, we have understood the matching of Wilson loop eigenvalue densities in both 
the strong-coupling and weak-coupling limits, and shown it is violated away from these limits. 
These violations imply that the Gross-Witten transition will not be matched in 2+1 and 3+1 
dimensions. The Durhuus-Olesen non-analyticity will be reproduced for bare loops in 2+1 
and 3+1 dimensions, but at an ultraviolet length-scale, and it will not be reproduced for 
smeared loops. 
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A Matching of eigenvalue spectra in the strong— coupling 
expansion 

As shown in Section 12. H matching the traces of two Wilson loops according to eq. ([1]) leads 
to matching of the traces of all powers of the loops, and hence to matching of the spectra, to 
leading order in the strong-coupling expansion. In this Appendix we consider the next order, 
described by tilings with one "bump" , and show that the matching also works there. 

We consider first the plaquette in SU{2). The strong-couphng expansion for the trace in 
D = 1 + 1 dimensions is 

Tr(P) = aiA + a:,(3l + a^(3l + (9) 

and in = 2 + 1, 3 + 1 dimensions it is 

Tr(P) = a^l3D' + a^l3l, + ^13%, + 0{l3l,). (10) 

The /3 and (3^ terms are the same in both cases since they correspond to only planar tilings, 
but the 13^ terms are different since they can include a "bump" in 2+1 or 3+1 dimensions. 
Similarly, the strong-coupling expansion of the trace of the plaquette squared is 

Tr(p2) = c^(3l + c^(3t + c^Pl + 0{(3l) (11) 

in 1+1 dimensions, and 

Tr(P^) = C2l3l, + C4/5^, + d,l3% + 0(/3« ,) (12) 
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in D' dimensions. The first term is rather than 0{(3), since we need to tile the doubly- 

winding loop with two plaquettes, the term has 4 plaquettes, the /3® term has either 6 
plaquettes in the plane or 1 in the plane and 5 in a bump, and so on. 

The matching condition ([T]) is that ([9]) and ( fTOl) are equal. We can expand 13 d' in the 
condition gives 

(3n'=f3, + '^^f3l + 0{f5l). (13) 
ai 

Substituting this into ( |T2l) we find the values of Tr(P^) will match in 1+1 and D' dimensions 
if the coefficients in the expansion are related as follows: 

4 = 06 + 2^(65-05). (14) 
ai 

The first term on the right-hand side is the 0{j3i) term in ( |TT|) . which is the same as the 
contribution from fiat tilings to d^. The second term should then give the contribution from 
tilings with a bump if (IMI) is to be satisfied. Now 65 — 05 is the difference between the C(/3^) 
contributions to (fTOj) and ([9]) i.e. it is the contribution from a "bump" in a single sheet. The 
factor ^ accounts for the different multiplicative factors we require when tiling a doubly- 
winding loop with two sheets rather than a singly-winding loop with one sheet. Finally the 
factor of 2 accounts for the fact that the "bump" can occur in either of the two sheets. Putting 
all these together gives exactly the correct expression for d^, so f|T^ is indeed satisfied. 

Generalising this to the trace of the /c'th power of the Wilson loop Tr(P'^) is straightfor- 
ward. It is also straightforward to generalise to Wilson loops of area A. Finally, changing 
N changes the number of plaquettes we can add to a tiling at one site from 2 to A^; this 
will change the 0{[3^) terms in (ED and ^ to 0(/?^+^), but it will not affect the "bump" 
contribution at 0{[3^) since this is singly-tiled everywhere. Hence the proof above will also 
work for all values of A^. 

So the first term in the strong-coupling expansion corresponding to tilings in which one 
of the sheets has a "bump" will indeed give matching of traces of powers of the Wilson loop 
if the matching condition ([T]) is satisfied. Hence the eigenvalue densities will also match up to 
that order in the strong-coupling expansion. 
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Plaquette in SU{2) 





Lattice size 






k' 


Fk' 


0.07 


42 
43 
44 


0.281 

0.2807 

0.2808 


0.28095 

0.28078 
0.28082 


2 


-0.495066(20) 
-0.495087(13) 
-0.495087(8) 


0.3 


42 
43 

1^ 


1.279 

1.2605 

1.212 


1.27893 
1.26049 

1.21110 


2 


-0.407127(32) 
-0.407189(37) 

-0.107195(15) 


0.4 


42 
43 
44 


1.803 
1.729 
1.648 


1.80340 
1.72878 
1.64803 


3 


-0.351050(35) 
-0.351121(18) 
-0.351256(25) 


0.5 


42 
43 

44 


2.444 
2.235 
1.996 


2.44677 
2.23511 

1.99071 


3 


-0.39550(17) 

-0.39629(7) 

-0.39769(24) 


0.6 


8^ 
83 
8^ 


3.31 
2.83 
2.29 


3.3114 
2.8279 
2.2933 


3 


-0.39758(18) 
-0.40010(14) 
-0.40638(13) 


0.7 


8^ 

83 
8^ 


4.64 
3.67 
2.804 


4.6440 
3.6734 

2.8033 


4 


-0.40855(14) 

-0.41185(6) 

-0.41736(9) 


0.8 


82 
83 
8^ 


7.17 

5.3 

4.0 


7.1940 
5.3006 
4.0062 


5 


-0.41510(35) 
-0.41921(14) 
-0.42260(20) 


0.9 


82 
83 
8^ 


14.7 
10.27 
7.73 


14.719 

10.253 
7.729 


7 


-0.43006(34) 
-0.43298(17) 
-0.43448(15) 


0.97 


82 
83 

8^ 


49.4 
33.5 
25.2 


49.598 
33.545 

25.218 


12 


-0.44220(14) 
-0.44332(5) 

-0.11385(5) 


0.99 


122 
123 
12^ 


149.01 
100.15 
75.2 


149.534 
100.211 
75.194 


21 


-0.444903(84) 
-0.445286(19) 
-0.445457(33) 



Table 1: Fourier components Fk' at the minimum k' of the Fourier spectrum of the eigenvalue 
distribution, for plaquettes in SU{2). 
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2x2 Wilson loops in SU{2) 



{Uw) 


Lattice size 






k' 




0.04 


8^ 


2.0883 


2.08746 


2 


-0.499370(10) 




83 


1.9434 


1.94312 




-0.499358(13) 




8^ 


1.7915 


1.79090 




-0.499344(10) 


0.07 


82 


2.555 


2.55399 


2 


-0.497887(13) 




83 


2.273 


2.27472 




-0.497781(18) 




8^ 


1.985 


1.98421 




-0.497697(25) 


0.15 


12^ 


3.552 


3.55501 


2 


-0.487888(29) 




123 


2.877 


2.87655 




-0.487604(24) 




12^ 


2.233 


2.23233 




-0.487568(34) 


0.3 


122 


5.438 


5.4385 


2 


-0.435366(23) 




123 


3.85 


3.8501 




-0.435169(22) 




12^ 


2.646 


2.6459 




-0.436626(77) 


0.4 


122 


7.019 


7.0179 


3 


-0.376907(63) 




123 


4.642 


4.6417 




-0.377007(39) 




12^ 


3.119 


3.1167 




-0.377989(79) 


0.5 


122 


9.13 


9.1355 


3 


-0.43410(7) 




123 


5.7 


5.7104 




-0.43470(8) 




12^ 


3.8 


3.8015 




-0.43584(19) 


0.6 


122 


12.22 


12.2357 


3 


-0.44130(13) 




123 


7.285 


7.2886 




-0.44210(8) 




12^ 


4.82 


4.8230 




-0.44279(5) 


0.7 


122 


17.3 


17.3138 


4 


-0.43348(21) 




123 


9.87 


9.8833 




-0.43417(14) 




12^ 


6.54 


6.5368 




-0.43445(10) 


0.8 


122 


27.39 


27.3759 


5 


-0.432074(48) 




123 


15.08 


15.0791 




-0.432562(47) 




12^ 


9.97 


9.9604 




-0.432752(37) 


0.9 


16^ 


57.36 


57.4346 


7 


-0.439275(47) 




163 


30.66 


30.6838 




-0.439529(36) 




16^ 


20.23 


20.2416 




-0.439604(42) 


0.97 


242 


197.5 


197.411 


12 


-0.445541(25) 




243 


103.5 


103.483 




-0.445624(14) 



Table 2: Fourier components Fk' at the minimum k' of the Fourier spectrum of the eigenvalue 
distribution, for 2 x 2 Wilson loops in SU{2). 
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4x4 Wilson loops in SU (2) 



{Uu,) 


Lattice size 




/^rcw 


k' 


Fki 


0.04 


162 


7.933 


7.93321 


2 


-0.499699(10) 




16^ 


4.497 


4.49635 




-0.499642(12) 




16^ 


2.636 


2.63619 




-0.499688(9) 


0.07 


162 


9.506 


9.50617 


2 


-0.498705(8) 




16^ 


5.067 


5.06701 




-0.498632(10) 




16^ 


2.8908 


2.89075 




-0.498708(9) 


0.2 


162 


15.4 


15.4000 


2 


-0.479298(11) 




16^ 


7.1843 


7.18433 




-0.479154(8) 




16^ 


3.9905 


3.99028 




-0.479322(11) 


0.4 


162 


26.69 


26.6860 


3 


-0.379362(15) 




16^ 


11.219 


11.2183 




-0.379351(18) 




16^ 


6.231 


6.23144 




-0.379401(20) 


0.5 


162 


35.124 


35.1209 


3 


-0.437299(16) 




16^ 


14.224 


14.2245 




-0.437275(17) 




24^ 


7.9327 


7.93339 




-0.437364(12) 


0.7 


242 


67.74 


67.7909 


4 


-0.435227(19) 




243 


25.974 


26.0474 




-0.435269(21) 




24^ 


14.563 


14.5656 




-0.435259(16) 



Table 3: Fourier components F^/ at the minimum k' of the Fourier spectrum of the eigenvalue 
distribution, for 4 x 4 Wilson loops in SU{2). 




5 10 15 20 

k 



Figure 1: Fourier components Fk for the plaquette in SU{2) with {up) = 0.5 on a 42 lattice. 
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Plaquetl 


e in SU{3) 








Lattice size 




(3tcw 


k' 




0.07 


42 


1.15 


1.15022 


2 


-0.062686(18) 




43 


1.15 


1.15038 




-0.062663(25) 




44 


1.152 


1.15112 




-0.062670(20) 


0.3 


42 


4.268 


4.2690 


2 


-0.165392(13) 




43 


4.187 


4.1864 




-0.165408(13) 




44 


4.102 


4.1017 




-0.165438(13) 


0.5 


42 


7.3 


7.2996 


2 


-0.11133(6) 




43 


6.484 


6.4920 




-0.11246(9) 




44 


5.51 


5.4985 




-0.11543(27) 


0.7 


82 


12.94 


12.9408 


3 


-0.16766(8) 




8=^ 


9.91 


9.9328 




-0.17274(6) 




8^ 


7.54 


7.5356 




-0.17800(8) 


0.97 


8^ 


132.8 


132.798 


11 


-0.193903(31) 




8=^ 


89.6 


89.597 




-0.194490(31) 




8^ 


67.4 


67.338 




-0.194797(34) 



Table 4: Fourier components F^i at the minimum k' of the Fourier spectrum of the eigenvalue 
distribution, for plaquettes in SU{3). 



Plaquette in SU{6) 





Lattice size 






k' 


Fk' 


0.3 


42 


20.78 


20.7839 


2 


-0.031469(28) 




43 


20.42 


20.4343 




-0.031658(51) 




44 


20.09 


20.0838 




-0.031747(26) 


0.435 


42 


29.7 


29.355 


2 


-0.04909(4) 




1^ 


21.15 


21.151 




-0.05650(12) 


0.5 


42 


34.0 


33.970 


2 


-0.03759(8) 




43 


29.7 


29.882 




-0.04316(6) 


0.52 


42 


35.5 


35.575 


2 


-0.03020(8) 




44 


24.15 


24.129 




-0.05103(24) 


0.7 


8^ 


58.0 


57.920 


3 


-0.09947(5) 




8^ 


43.8 


43.788 




-0.10963(5) 




8^ 


33.1 


33.160 




-0.11784(9) 


0.97 


82 


582.3 


582.73 


10 


-0.13828(11) 




8=^ 


391.0 


391.28 




-0.13952(7) 




8^ 


294.2 


294.59 




-0.14012(7) 



Table 5: Fourier components F^' at the first minimum k' of the Fourier spectrum of the 
eigenvalue distribution, for plaquettes in SU{6). 
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2x2 Wilson loops in SU{3) 



(v ) 


T li^f f i pp til 7P 

J-JCXIj OIZjC 


B 

yrun 


B 


k' 




yj. yjrt 


8^ 




fi 3875 


2 


_n 03891 




83 
o 


o.oo 


^ ^941 

tJ. OZi'rt J. 




u.uooo MM 




8^ 




^ 91 01 






0.15 


12^ 


10.12 


10.1213 


2 


-0 130071f21) 




123 


7.857 


7.8574 




-0.129931(46) 




12^ 


5.774 


5.7742 




-0.130398(56) 


0.5 


122 


24.86 


24.8285 


2 


-0.14448(8) 




123 


15.17 


15.1684 




-0.14497(13) 




12^ 


10.1 


10.1042 




-0.14548(13) 


0.9 


122 


153.5 


153.412 


6 


-0.194108(25) 




123 


81.7 


81.767 




-0.194218(56) 




12^ 


53.97 


53.959 




-0.194307(62) 



Tabic 6: Fourier components Fk' at the minimum k' of the Fourier spectrum of the eigenvalue 
distribution, for 2 x 2 Wilson loops in SU{3). 



2x2 Wilson loops in SU{6) 



{Uw) 


Lattice size 






k' 


Fk' 


0.04 


82 


30.18 


30.1791 


2 


-0.008415(31) 




83 


27.47 


27.4617 




-0.008467(34) 




8^ 


24.42 


24.4228 




-0.008630(52) 


0.15 


122 


45.8 


45.8021 


2 


-0.066150(48) 




123 


34.8 


34.8034 




-0.067323(34) 




12^ 


25.035 


25.0292 




-0.069585(53) 


0.5 


122 


109.5 


109.584 


2 


-0.10342(6) 




123 


66.4 


66.259 




-0.10453(15) 




12^ 


44.25 


44.2464 




-0.10534(13) 


0.9 


122 


672.4 


672.859 


6 


-0.137916(64) 




123 


358.0 


357.767 




-0.138330(41) 




12^ 


236.0 


236.098 




-0.138396(50) 



Table 7: Fourier components F^' at the first minimum k' of the Fourier spectrum of the 
eigenvalue distribution, for 2 x 2 Wilson loops in SU{6). 
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1 1 r 




-3-2-10123 

a 

Figure 2: Eigenvalue densities for the plaquette in SU{2) with {up) = 0.7 in D = 1 -|- 1 (sohd 
hne), D = 2 + 1 (dashed) and D ^ 3 + 1 (dotted). 
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Figure 4: Difference between F^/ in 1+1 dimensions and D' dimensions, for = 2 + 1 
and D' — 3 + 1 (x), as a function of the trace (up). All for the plaquette in SU{2). 
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Figure 5: Difference between Fk' in 1+1 dimensions and D' dimensions, for D' = 2 + 1 (+) 
and D' = 3 + 1 (x), as a function of the trace All for the 2x2 Wilson loop in SU{2). 
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Figure 6: Difference between F^/ in 1-1-1 dimensions and D' dimensions, for = 2 -|- 1 (-|- 
and D' — 3 + 1 (x), asa function of the trace (uu,)- All for the 4x4 Wilson loop in SU{2). 
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Figure 7: Difference between Fk' in 1+1 dimensions and D' dimensions, for = 2 -|- 1 (-|-) 
and = 3-1-1 (x), as a function of the trace {up). All for the plaquette in SU{3). 
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Figure 8: Difference between F^/ in 1+1 dimensions and D' dimensions, for = 2 + 1 
and D' — 3 + 1 (x), as a function of the trace (up). All for the plaquette in SU{6). 
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Figure 9: Difference between Fk' in 1+1 dimensions and D' dimensions, for D' = 2 + 1 (+) 
and D' = 3 + 1 (x), as a function of the trace All for the 2x2 Wilson loop in SU{3). 
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Figure 10: Difference between F^/ in 1+1 dimensions and D' dimensions, for D' = 2 + 1 (+) 
and D' — 3 + 1 (x), asa function of the trace (uw). All for the 2x2 Wilson loop in SU{6). 
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